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I.  Introduction  and  Objectives 

» 

The  purpose  of  this  report  is  to  present  work  performed  under  this  AFOSR  contract. 
As  has  been  repeatedly  discussed  elsewhere  (Ref.  1-5)  the  thrust  of  this  project 
is  triple  pronged:  1)  to  calculate  on  a  consistent  and  comprehensive  basis 
t  basic  energy  quantities  for  the  triatomic  molecules  H20,  C02,  N02,  N20, 

and  O3  and  their  positive  and  negative  ions.  These  quantities  include  the 
vertical  spectra,  the  ionization  potentials  and  the  electron  affinities 
of  these  systems.  As  part  of  this  effort  portions  of  potential  energy  surfaces 
>  are  computed  and  used  to  predict  vibrational  spectra  and  to  interpret  the 

mechanisms  of  chemical  reactions  of  interest.  2)  second  the  development 
of  two  new  methods  for  treating  electronic  structures  will  be  discussed. 

3)  the  third  topical  area  of  this  research  involves  the  study  of  specific 
'  atmospheric  problems  involving  these  molecular  systems  as  they  are  brought 

to  our  attention  through  the  literature  or  directly  by  experimentalists. 

In  this  section  a  molecule  by  molecule  summary  will  be  given. 

This  report  will  follow  the  above  structure.  First,  overall  status  will 
be  given,  then  progress  on  two  new  method  development  projects  (MCSCF  and 
Differential  Correlation  Energy  (DCE)  will  be  summarized,  next  a  molecule 
by  molecule  summary  summary  will  be  given. 
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II.  Overall  Status  of  Molecular  Project 

The  overall  status  of  the  triatomic  project  is  summarized  in  Table  I.  This 
is  a  cumulative  table  with  results  obtained  during  the  current  contract 
period  (March  1,  1978  -  May  30,  1980),  given  in  parenthesis. 

Reference  to  papers  published  covering  work  in  Table  I  are  given  in  Appendix 

I. 
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III.  The  Differential  Correlation.  Energy  Method 

a.  Introduction 

Although  potential  curves  of  chemical  accuracy  can  now  be  computed  routinely, 
the  calculation  of  electron  affinities  and  ionization  potentials  still  presents 
a  considerable  challenge.  Methods  which  have  been  used  to  calculate  electron 
affinities  and  ionization  potentials  are  extensive  configuration  interaction, 
equations  of  motion,  pair  energy  schemes,  Hess  cycles,  perturbation  theory, 

Koopmans  theorem  and  Hartree  Fock  methods. 

The  Hartree-Fock  model  (and  Koopman's  theorem)  usually  yield  electron  affinities 
and  ionization  potentials  which  are  in  error  by  several  ev,  however  through  fortuitous 
cancelation  of  basis  set  and  correlation  errors,  good  agreement  with  experiement 
is  sometimes  obtained.  The  errors  in  other  methods  are  usually  smaller  however 
only  extensive  configruation  interaction  or  Hess  cycle  calculations  have  yielded 
quantitative  results. 

In  this  section  we  present  a  configuration  selection  scheme  designed  to  only  in¬ 
clude  those  excitations  which  directly  contribute  to  the  correlation  energy 
difference  between  the  neutral  and  ionic  system.  The  orbitals  for  each  con¬ 
figuration  are  evaluated  via  the  MCSCF  method.  We  call  this  scheme  the  Differential 
Correlation  Energy  (DCE)  method. 

b.  The  Differential  Correlation  Energy  Method 

We  are  currently  developing  a  new  method  for  choosing  configurations  in  a 
MCSCF/CI  calculation  to  give  accurate  (within  0.1  eV)  electron  affinities  and 
ionization  potentials  using  the  smallest  number  of  configurations  possible. 

The  key  to  this  method  is  the  determination  of  the  configurations  which  con¬ 
tribute  to  the  differential  correlation  energy  (DCE)  associated  with  the 
addition  or  removal  of  an  electron.  It  is  our  aim  to  develop  a  prescription 
for  the  choice  DCE  configurations  analogous  to  the  OVC  prescription  6  for  the 
choice  of  configurations  to  give  accurate  (+  .1  eV)  bond  dissociation  energies. 

THE  OVC  method  concentrates  only  on  the  changes  in  correlation  energy  associated 
with  bond  formation. 

The  DCE  method  takes  advantage  of  some  simple  properties  of  correlation  energy. 


(9) 


1)  Correlation  energy  is  very  sensitive  to  electron  occupations,  i.e.,  new 
pairs,  orbital  occupations. 

2)  Correlation  energy  is  relatively  insensitive  to  details  of  orbital  shapes. 

3)  Pair  energy  schemes  work. 

In  light  of  the  above. 

1)  It  is  reasonable  to  evaluate  only  the  correlation  energy  difference 
between  neutrals  and  ions  or  different  states. 

2)  Major  di fferences  come  from  new  pairs  formed  and  cross  shell  correlation 
of  extra  electrons  with  remaining  ones.  (Active  orbital  contribution). 

3)  Minor  differences  some  from  pair  correlation  energy  contribution  due  to 
availability  of  open  orbital.  (Excluded  space  effect.) 


c.  Formal  Development  for  Electron  Affinities  and  Ionization  Potentials 

Below  we  give  specifically  the  formulation  of  ideas  proposed  for  the  calculation 
of  electron  affinities  by  the  DCE  method. 


In  order  to  discuss  electron  affinity  or  ionization  potential,  it  is  convenient 
to  use  a  second-quantized  treatment.  Thus  we  introduce  the  usual  electron  annihilation 
and  creation  operators  ai  and  aj,  respectively.  The  index  i  serves  to  distinguish 
among  different  spin-orbitals  in  our  ordered  base.  The  operators  satisfy  the 
anticommutation  relations 

{ai ,aj )  =  taj’  aj>  =  0  ’ 

■  5ij  • 

Now  if  <J>  denotes  an  electronic  state  function  (in  Fock  space),  a*  $  will 
describe  the  same  state  function  with  an  additional  electron  in  the  i-th  orbital. 

This  assumes  that  the  i-th  orbital  is  unoccupied  in  <&.  For  simplicity,  we  shall 
assume  that  <J>,  and  consequently  aj$,  represent  single-configuration  state  functions. 
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The  Hamiltonian  can  be  written 

H  =  Eij  al  hu  aj  +  *  Eij  EW  WiJ,EVj 

This  can  be  partitioned  according  to 

H  =  H  +  H 

0 


where 


Ho  ■  vt  hh  aj  • 


H  "  Eijai  ^ij'hij^  aj  +  15  Ei,j  Ek£aiak9i;j,k£.a£aj  . 

We  have  not  specified  our  orbital  base  in  detail,  so  we  are  free  to  assume 


that  it  is  chosen  so  that 


hij  Sij  £i  ■ 

Then  4>  and  a|$  are  eigenvectors  of  Hq;  that  is, 
Ho  «  -  E  *  , 


Vi4,  =  !Eo  +  Ei>  4  • 


where 


I  •  .  £  • 
ie$  l 


We  have  some  freedom  in  the  choice  of  h°;  for  example,  we  can  choose  h°  so 
as  to  give  Hartree-Fock  orbitals  (for  the  state  $)  as  its  eigenfunctions. 
Note  that  in  this  case  Eq  is  not  the  Hartree-Fock  energy  for  the  state; 
the  Hartree-Fock  energy  is  given  by  <$]H|$>. 


No  matter  what  the  choice  of  h°  we  use,  4>  will  not  be  an  eigenvector  of  H 
because  of  electron  correlation  effects,  and  the  same  observation  applies 
also  to  aj<i>.  We  can  compute  these  correlation  effects  approximately  by 
application  of  Rayleigh-Schrodinger  perturbation  theory,  with  Hq  as  the 
unperturbed  Hamiltonian  and  H1  as  the  perturbation. 


(ID 


This  yields  to  first  order. 


f  5  $>  +  R$  , 

V  =  +  R  at<j)  , 


where 


R  *  £io,ivj  4  aj/(Ej-Ei> 

+wijEk£,  c1+e^Sj+S£  aK9ij,wataj/(£i+£j'crEk!  • 


To  second  order  the  energies  are  given  by 
E  =  <4>  |  H 1 4»  +  <$1 H 1  R  |4»  , 

Ee  =  <<t>  { a  ^  Ha^  ( +  <<?  ( a  ^  H '  Rajf'Jo  . 

Since  we  have  a.a+  ®  =  4>,  the  electron  affinity  is  given  to  second  order 

i  7 


Ee-E  =  e.+<$|  [a^H'J  a*|$> 
+«t>|[a.,H'J  Ra*  +  H'  [a^Rj  a{|$> 


From  this  expression  we  see  that  the  computed  electron  affinity  remains 
unaltered  to  this  order  if  we  modify  H’  by  adding  to  it  any  Hermitian 
operator  that  commutes  with  a^.  Hence  we  can  us<*  in  place  of  H'  the 
expression 


H"  =  a.  H+  +  HJi. 


where 


Hi  - l,  [hij-hij+  zk  4  at9ij;k*]aj' 

H-’Vjhr  h5(  *  £wakVji,w]  . 


(12) 


In  other  words,  any  part  of  H"  that  does  not  reference  the  orbital  i  does 
not  contribute  to  the  computation  of  the  electron  affinity. 

It  is  not  advisable  to  take  the  above  argument  too  literally.  For  example, 
it  does  not  indicate  that  many-body  perturbation  theory  will  give  a  good 
result  for  the  electron  affinity  if  carried  to  second  order  only.  Higher 
order  effects  arelikely  to  prove  important.  Many  of  these  effects  are 
automatically  included  in  a  SCF  or  MCSCF  treatment,  so  we  shall  not  pursue 
such  questions  further.  What  remains  is  a  guide  in  attempting  to  answer 
the  question:  What  configurations  are  needed  to  handle  electron  correlation 
effects  that  directly  affect  the  calculation  of  electron  affinity?  From  the 
above  discussion,  the  configuration  lists 

$,  aja^a^  $  ,  any  j,  k,  £ ,  f,  i, 

a*  4>,  ,  any  j ,  k,  l ,  f,  i , 

appear  as  the  answer  to  this  question.  In  essence,  this  tells  us  that 
for  the  state  without  the  added  electron  we  should  include  the  dominant 
(i.e.,  Hartree-Fock)  configuration  and  all  single  and  double  replacements 
which  contain  the  orbital  of  the  electron  to  be  added.  For  double  replace¬ 
ments,  we  note  that  in  addition  to  the  spin  orbital  i,  the  configuration 
may  also  contain  another  orbital  from  this  shell.  For  the  state  with  the 
added  electron,  we  should  again  include  the  dominant  configuration  and  any 
single  and  double  replacements  that  are  generated  by  el iminatinq  the  spin 
orbital  of  the  added  electron  in  the  dominant  configuration. 
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Some  obvious  refinements  can  be  added  to  this  scheme:  single  replacements 
from  the  dominant  configuration  are  unlikely  to  be  important,  since 
Brillouin's  theorem  applies  to.  an  SCF  approach,  and  replacements  which 
involve  excitation  of  a  core  orbital  are  probably  negligible  because  a 
large  excitation  energy  is  involved  (the  core  orbitals  can  be  defined  as 
those  orbitals  to  which  such  an  argument  applies).  We  also  must  be 
concerned  with  the  need  for  a  multiconfiguration  reference  state. 

A  similar  analysis  is  being  performed  for  the  calculation  of  vertical 
spectra. 

d.  Appl ications 

Using  the  DCE  configuration  selection  rules  MCSCF  calculations  have  been 
performed  on  CN->£N~  and  C02-*C02. 

2 

For  CN  we  must  include  all  correlation  within  the  5o  orbital  and  all  cross 

2  2  2 

shell  correlation  of  5o  with  3o  ,  4o  ,  and  li.  Using  3  new  MCSCF  correlating 
orbitals,  6o,  7a,  2tt  this  yields  25  configurations. 

For  CN  the  DCE  perscription  tells  us  to  include  all  cross  shell 
correlation  involving  the  5a  orbital  and  all  excitations  into  the 
half  filled  5a  orbital.  This  yields  41  configurations.  Results 
of  these  calculations  are  presented  in  Table  II. 

For  C02  we  must  include  all  valence  excitations  into  the  6a^  orbital,  this 
yields  67  configurations.  For  CO^  we  must  include  all  cross  shell  correlation 
of  6aielectron  with  valence  shell  and  valence  excitations  into  68^  orbital. 
Results  on  CO^CQ^  are  presented  in  Table  III. 

These  results  indicate  that  the  DCE  rules  may  provide  a  reasonable  starting 
point  for  the  MCSCF  procedure  which  can  then  be  followed  by  a  larger  Cl  for 
for  truly  quantitative  results. 
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TABLE  II 

DCE  METHOD  APPLIED  TO  CN  CN" 

MCSCF/CI 


State/Enerqy 

15/26 

Single/Conf . 

SCF 

AE 

CN"  R=2. 117 

-92.3724551 

-92.321633 

-92.3221717 

.0508 

"  2.217 

-92.3812117 

-92.326182 

-92.3267819 

.0550 

"  2.317 

-92.3780383 

-92.318317 

-92.3189866 

.0597 

CN( 19MC )2 . 117 

not  converged 

-92.2045157 

"  2.217 

I*  II 

-92.2027169 

"  2.317 

-92.2510795 

-92.172356 

-92.1890972 

.0787 

CN(fromCN") 

2.117 

-91.9523405 

-91.905750 

-92.2045157 

.0466 

2.217 

-91.9451404 

-91.897083 

-92.2027169 

.0491 

2.317 

-91.9445647 

-91.895993 

-92.2090972 

.0486 

CN  (17MC)2. 117 

-92.2317141 

-92.203380 

-92.203380 

.0312 

2.217 

-92.2402757 

-92.201029 

-92.2027169 

.0392 

2.317 

-92.2358326 

-92.186712 

-92.1890972 

.0491 
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TABLE  III 


DCE 

METHOD 

APPLICATION  TO  C0^C02 
(0oco  *  135.77°,  Rco  =  2.3734  bohrs) 

Configurations  for  C02 
Orbitals 

1  a i  1 b2  2ai  3a i  4b2  4aj  2b2  5ai  3b2  6ai  7ai  5b2  1 b i  la2  2bi 
222222222100220 

1)  Include  all  cross  shell  correlation  of  6ai  electron  with  valence  electrons 

2)  Include  all  valence  excitations  into  6a i  orbital. 

E  =  -187.56976  a.u.  19MCSCF/70CI 

Configurations  for  C02 

222222222000220 
1)  Include  all  valence  excitations  into  6a t  orbital 
E  =  -187.52993  a.u.  19MCSCF/67CI 

AE  =  1 .08ev 

AE  =  .3ev  (basis  set  error  in  C02) 

Vertical  Ionization  Potential  of  C02  -  1.4ev 
Experimental*  IP  of  COl  1.4+. 3ev 
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IV.  A  Multidimensional  Newton-Raphson  MCSCF  Method 

The  merits  of  the  MCSCF  approach  to  the  calculation  of  electronic  wave 
functions  for  atoms  or  molecules  are  well  documented^and,  in  fact,  the 
method  has  become  the  standard  starting  point  for  most  researchers.  The 
principal  limitation  of  this  approach  has  been  the  frequent  difficulty  in 
obtaining  converged  solutions  for  the  MCSCF  wave  function.  A  new  MCSCF 
method 'under  development  by  J.  Detrich  and  C.C.J.  Roothaan  is  expected 
to  eliminate  this  difficulty. 

This  method  applies  whenever  the  spin-orbitals  (SO)  are  put  forward  as 
finite  expansions  in  terms  of  a  basis  set.  In  the  MCSCF  approach  one  also 
constructs  term  state  functions  (TSF)  as  linear  combinations  of  configuration 
state  functions  (CSF);  the  latter  are  defined  as  symmetry-adapted  linear 
combinations  of  Slater  determinants  arising  from  a  single  configuration  of 
orbital  assignments.  A  complete  set  of  trial  SO's  can  be  put  forward,  which 
spans  the  same  linear  space  as  the  one-electron  basis  functions,  and, 
analogously  a  set  of  trial  TSF's  spanning  the  same  linear  space  as  the  set  of 
CSF's  can  be  developed. 

The  method  deals  in  terms  of  the  unitary  transformations  of  the  reference 
SO’s  and  TSF's  which  will  give  optimized  SO's  and  TSF's.  These  unitary 
transformations  are  expressed  as  exponentials  of  anti-Hermitian  operators. 

In  terms  of  the  reference  SO's  and  TSF's,  the  variational  energy  expression 
is  of  the  form 

where  the  x.'s  are  the  independent  matrix  elements  of  the  anti-Hermitian 
operators.  The  quantities  E,  E.  and  E. .  are  evaluated  by  rather  elaborate 

I  J 

contractions  involving  integrals  over  basis  functions,  occupation  numbers 
and  coupling  constants,  and  the  expansion  coefficients  representing  the  set  of 
reference  SO's  and  TSF's. 

Approximating  this  expression  by  dropping  all  terms  of  third  and  higher 
order,  the  extremum  of  E'  with  respect  to  all  the  variables  x..  occurs  when 


they  satisfy 


■  -E>- 


Once  the  solution  of  this  equation  has  been  obtained,  new  reference  SO's  and 
TSF's  are  easily  computed,  thus  completing  an  iterative  MCSCF  cycle. 
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This  method  is  seen  to  be  a  multidimensional  generalization  of  the  Newton- 
Raphson  iterative  process,  and  it  converges  quadratically  in  the  general 
case.  The  availability  of  such  a  method  will  substantially  enhance  the 
range  of  problems  where  the  advantages  of  the  MCSCF  approach  can  be 
exploited.  Even  in  cases  where  other  MCSCF  methods  suffice,  the  availability 
of  a  method  that  virtually  guarantees  convergence  in  a  few  iterations  is 
expected  to  significantly  reduce  the  man-hours  required  to  perform  a  given 
series  of  calculations.  These  advantages  indicate  that  the  new  method  will 
become  the  standard  vehicle  for  molecular  orbital  calculations  in  the  basis- 
set  expansion  regime. 

A  preliminary  implementation  of  this  method  has  been  coded  for  atoms  under 
the  auspices  of  the  National  Resource  for  Computation  in  Chemistry  by 
Prof.  Roothaan.  A  comprehensive  description  of  the  formalism  and 

analysis  for  the  method  is  currently  in  preparation. 


Professor  C.  C.  J.  Roothaan  is  a  technical  consultant  on  self  consistent 
field  techniques  to  SAI. 
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Molecule  by  Molecule  Summar; 


CO2  and  its  Ions 

Work  on  C02  and  its  ions  during  the  contract  period  consisted  primarily 
of  trial  calculations  intended  to  help  establish  and  efficient  and  reliable 
approach  to  the  calculation  of  electron  affinities  and  ionization  potentials. 
As  discussed  in  Section  III  above,  this  goal  has  motivated  development 
of  the  Differential  Correlation  Energy  (DCE)  method.  The  calculations 
for  C02  and  C02  summarized  here  serve  as  a  test  for  the  method. 


For  both  C02  and  C02  the  (3s, 2p)  contracted  Cartesian  Gaussian  basis  of 

g 

Dunning  was  adopted.  This  choice  permits  trial  calculations  to  be  economically 

performed,  but  involves  some  sacrifice  in  accuracy  of  the  orbital  description. 

As  an  aid  in  assessing  such  effects,  SCF  vertical  spectra  for  C02  were  compu;.;'.' 

with  the  (3s,  2p)  basis  and  compared  with  previous  results  using  a  more  flexible 
.  9 

basis  at  the  experimental  linear  equilibrium  geometry  (Table  IV)  and  a  bent 

geometry  (Table  V).  In  general,  the  SCF  vertical  spectra  for  C02  are  in  good 

agreement,  the  only  exception  being  in  the  case  of  the  V  state,  which  is 

predominantly  Rydberg  rather  than  valence  in  character,^  and  hence  not  well 

described  by  a  basis  set  which  does  not  contain  diffuse  functions.  This  indicates 

that  the  correlation  effects  of  interest  are  reasonably  well  described  in  terms  of 

the  (3s,  2p)  basis.  Another  indication  of  the  utility  of  the  (3s,  2p)  basis  comes 

from  the  results  of  an  SCF  geometry  scan  for^the  ;gr&und  state  of  C0£.  The 

equilibrium  geometry  is  found  to  be  0OCO  =  136°,  Rcq  =  2.37  bohr,  in’iood 

aggreement  with  the  determination  eoco  =  135°,  RcQ  =  2.35  bohr  by  Pacansky,  Wahlgren 
•  1  \ 

and  Bagus*  using  a  considerably  larger  basis. 


The  DCE  method  was  applied  to  compute  the  vertical  ionization  potential 
of  C02  at  the  SCF  equilibrium  geometry  60C0  =  135.77°,  Rco  =  2.3734  bohrs. 
In  addition  to  the  orbitals  occupied  in  the  dominant  ground  state  C02  con¬ 
figuration,  three  correlation  orbitals,  7a155b2,  2bp  were  introduced. 

The  list  of  configurations  developed  according  to  the  DCE  method  contained 


70  for  CO^  and  67  for  CO,,.  The  orbitals  were  taken  from  an  MCSCF  calculation 
for  COj  involving  the  19  most  significant  configurations  from  the  list 
of  70  OCE  configurations  for  COg.  Note  that  the  orbitals  from  CO2  calculation 
were  used  in  the  CO2  as  well;  this  was  the  only  obvious  way  to  preserve 
the  character  of  the  6a^  orbital  in  the  COg  calculation. 

The  energy  computed  from  the  67-configuration  interaction  calculation  for 
CO^  was  E  =  187.52993  hartrees;  from  the  70-configuration  interaction 
calculation  for  CO^  the  energy  E  =  -187.56976  was  obtained.  This  implies 
a  vertical  ionization  potential  for  CQ^  &E  =  1.G8  eV. 

As  noted  above,  the  basis  set  adopted  affects  the  accuracy  of  this  result. 
Since  the  effect  is  not  expected  to  depend  strongly  on  correlation,  one 
can  obtain  a  quantitative  estimate  by  comparing  the  SCF  result  for  the 
CO2  adiabatic  electron  affinity  computed  with  the  (3s, 2p)  basis,  -.65  eV, 
with  the  corresponding  SCF  result  from  Pacansky,  et.al .l1  -.38  eV.  This 
indicates  a  basis  set  error  of  about  .3  eV.  Revising  our  computed  result 
accordingly  gives  1.4  eV  for  the  vertical  ionization  potential  of  CO^. 

Strictly  speaking,  there  is  no  corresponding  experimental  determination, 

but  a  close  analogue  can  be  constructed  as  follows.  When  C02  is  distorted 

from  its  experimental  equilibrium  geometry  to  the  theoretical  CO^  equilibrium 

geometry  adopted  in  these  calculations,  the  CO2  potential  constructed  by 
12 

Suzuki  from  spectroscopic  data  gives  a  distortion  energy  of  2.0  eV. 

Adding  this  to  the  experimentally  derived  CO ^  electron  affinity  of  -  .6+  .2  eV 
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reported  by  Compton,  Reinhardt  and  Cooper  yields  1.4±  0.3  eV  for  the  C02 
vertical  ionization  potential. 

This  result  is  in  good  agreement  with  the  result  based  on  the  DCE  method. 


More  refined  computations  are’  necessary  to  confirm  our  estimate  of  basi 
set  error,  but  our  present  results  already  establish  the  DCE  method  as 
a  promising  new  approach. 


State 


TABLE  IV 

LINEAR  C02  VERTICAL  SPECTRUM 


r^=r2  =  2.1944  bohrs  6  =  180 


Configuration(s) 


Energy 


llTux  lirqx  2irux  lngy  XV  2*gy  — ^ 


4s, 3p  Energy 


xxLg+ 


2  0  -187.553005  0.00  0.00 


1  1 


1  1 


-187.292459  7.09  7.07 


-187.270395  7.69 


-187.270654  7.68 


-187.249106  8.27 


-187.249361  8.26  8.26 


-187.238979  8.54 


1  \  -187. 

0  ) 


279319  7.45 


■186.901777  17.72  14.15 


aCore  =  la  2  la  2  2o  2  3a  2  2a  2  4a  2  3a  2 
g  u  g  g  u  g  u 


1  a.u.  =  27.21070  eV 


'From  reference  9. 


TABLE  V 

BENT  C02  VERTICAL  SPECTRUM 
f^3s,2pj  Basis,  SCF 


rfr 

2~  2‘ 

1944 

bohrs 

8=130° 

State 

Configuration3 

Energy 

4s/3p  Energy' 

5d  i 

4b2 

6a1 

lbl 

la2 

2bl 

a.  u. 

eV 

eVb 

2 

2 

0 

2 

2 

0 

-187.462371 

0.00 

0.00 

2 

1 

1 

2 

2 

0 

-187.394377 

1.86 

2.00 

13A2 

2 

2 

1 

2 

1 

0 

-187.363691 

2.70 

2.76 

1XA2 

2 

2 

1 

2 

1 

0 

-187.352623 

3.00 

3.07 

11b2 

2 

1 

1 

2 

2 

0 

-187.  316653 

3.98 

4.05 

23B2 

2 

2 

0 

2 

1 

1 

-187.239182 

6.09 

6.13 

23A2 

2 

1 

0 

2 

2 

1 

-187.220412 

6.60 

6.70 

2:A2 

2 

1 

0 

2 

2 

1 

-187.208979 

6.91 

7.01 

2!B2 

2 

2 

0 

2 

1 

1 

-187.080537 

10.40 

10.39 

a 

Core= 

la^lb 

22a  2 
2  2al 

3a122b224a 

2  2 
iV 

b. 

la.  u. 

=  27. 

21070 

cReference  9 
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N:0+  x2n  MCSCF/CI  Calculation 

x  o 

A  geometry  scan  for  the  equilibrium  bond  lengths  of  the  linear  N20  X  it  state 
was  performed.  The  wave  functions  were  constructed  from  orbitals  expanded 
in  terms  of  a  basis  set  consisting  of  a  4s, 3p  Id  Gaussian  basis  on  each 
center.  The  MCSCF  calculations  consist  of  20  configurations,  and  the  sup¬ 
plementary  Cl  calculations  consist  of  99  configurations. 

The  results  for  the  linear  NNO  geometry  indicate  equilibrium  bond  lengths 
Rnn  =  2.126  a.u.,  R^q  =  2.367  a.u.  A  preliminary  fit  to  a  molecular  potential 
function  of  the  form 

»  ■  *  %  <4»nn>2  +  *  hio  <arno>2 

yields  =  19.563  md/A°,  =  5.144  md/A° 

A  more  elaborate  characterization  of  the  potential  surface  is  needed  in 
order  to  extract  spectroscopic  characteristics  such  as  vibrational  behavior. 

NOp"  Ground  State  Geometry  Scan 

A  preliminary  SCF-level  geometry  scan  of  the  N02“  ground  electronic  state 
has  been  performed.  They  include  the  normal  X*A^  C2v  state  of  equilibrium 
geometry,  the  peroxy  well,  and  the  energy-minimum  pathway  between  the  two 
geometries.  As  indicated  in  Table  VI,  the  geometries  of  the  C2v  and  peroxy 
Cs  minima  are  in  good  agreement  with  the  earlier  results  of  Pearson, 
et.  al.^  Hence  the  delineation  of  the  potential  surface  "between  these  two 
geometries  is  useful,  although  somewhat  crude.  It  would  appear  that  the  indi¬ 
cated  saddle  point  between  the  two  geometries  (Table  VI)  is  above  the  molecular 
dissociation  threshold,  so  that  isomerization  cannot  be  expected  to  take  place. 
Evidently  this  consistent  with  the  experimental  observations  that  bear  on  the 
existence  of  a  no2  n^roxy  well,  on  the  other  nano,  a  oreliminary  calculation 
cannot  be  relied  upon  with  regard  to  such  questions,  and  a  more  elaborate  set  of 

calculations  utilizing  multiple  configurations  and  a  larger  basis  set  is  planned. 
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TABLE  VI 


PRELIMINARY  SURVEY  OF  THE  N02  X1AJ  (V)  HYPERSURFACE 


SCF ,  (9s5p)/ 


a 


Critical 

ET  +  203 

V 

rno 

Roo 

aono 

aoon 

Geometry 

(a.u.j 

(eV) 

(A) 

(A) 

(deg.) 

(deg.) 

-1.0247 

0.0 

1.27 

2.17 

116.8 

31. 

(-1.0336) 

(0.0) 

(1.26) 

(117.0) 

S.P. 

-0.83 

5.3 

2.28 

2.06 

63.8 

82.6 

Ring 

(-0.8736) 

(4.4) 

(1.50) 

(58.5) 

Cs 

-0.9067 

3.2 

1.25 

1.51 

34. 

118.2 

(-0.9161) 

(3.2) 

(1.25) 

(1-49) 

(118.5) 

D(ON-O') 

4.0 

*The  values  in  parenthesis  are  (9s5p)/  5s3pJ 


results  from  reference  14. 
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The  wave  functions  utilized  in  this  set  of  calculations  are  based  on  the 

g 

(3s, 2p)  contracted  gaussian  basis  sets  of  Dunning.  In  support  of  these 
calculations,  the  l^Aj,  l^,  1^B2>  and  1^A2  states  of  NC^were 

computed  using  the  same  basis  set.  In  Table  VII  these  resultsare  compared 

*  Q 

with  earlier  results J  using  the  more  elaborate  4s, 3p  basis  set. 

3  1  3 

Similar  calculations  were  performed  for  the  lB^,  1  B2  and  1  states 

of  N02,  as  well  as  the  X^A^  ground  state.  These  are  of  some  interest  in 

their  own  right,  since  they  use  the  computed  ground  state  N02  equilibrium 

10 

geometry,  rather  than  the  N02  equilibrium  geometry  used  previously,  or 

15 

the  experimental  geometry  adopted  by  Benioff  .  Our  results  for  the  computed 

excitation  energies  are  given  in  Table  VUI.  Our  X1A1->11B2  excitation  energy 

appears  to  be  in  better  agreement  with  experiment  than  obtained  in  previous 
8  13 

work,  ’  but  more  elaborate  calculations  are  required  to  clarify  whether 
this  agreement  is  fortuitous. 
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TABLE  VII 


COMPUTED  SCF  ENERGIES  FOR  N02 


a 

State 

Configuration13 

Excitation 

3s2 £ 

Enerav  (eVV 
As3£ 

x  2a2 

4b2  ba} 

lb2  la2 

0.00 

0.00 

i  4a2 

4b2  ba^ 

lbj  lag  2b1 

2.05 

2.13 

1  4b2 

4b2  6ax 

lb2  la22b1 

2.10 

2.14 

12Bi 

4b2 

lb2  la22b1 

2.50 

2.50 

1  2A2 

4b2  6a2 

lb2  la2 

2.80 

2.83 

4b26a2 

2  2 
lb^la2 

3.82 

3.84 

a  All  states  at  bond  length  of  2.25  a.u.  and  bond  angle  of  134° 

k  The  omitted  orbitals  la^,lb2>  2a^,  3a^,  2b2,  4a^,  3b2,  5a^  are  doubly 
occupied  in  each  configuration. 

c  The  SCF  total  energy  of  the  X  2A,  state  is  -20 3.95597  a.u.  for  the  4s3p 
basis  set  and  -203.94646  a.u.  for  the  3s2p  basis  set.  The  difference  in 
eV  is  0.26. 
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.  TABLE  VIII 

EXCITATION  ENERGIES  FOR  N02 

b  Excitation  Energy (eV)  . 

State”  Configuration0  This  Workc  Experiment” 


4b2  6a^ 

lb? 

>•! 

0.0 

0.0 

4b2  6  a2 

lb? 

la2 

2b  l 

2.30 

- 

4b2  6  a2 

lb? 

la2 

2bl 

6.08 

~63 

4b2  6  a2 

Ibi 

2bi 

7.43 

- 

a)  All  states  at  bond  length  2.40  bohr  and  bond  angle  of  116.8° 

b)  The  omitted  orbitals  la1  lb2  3a^  2b2  4a1  3b,  5aj  are  doubly 
occupied  in  each  configuration. 

c)  Vertical 

d)  Adiabatic 

e)  W.  G.  Trawick  and  W.  H.  Eberhardt,  J.  Chem.  Phys.  22_,  1402  (1954); 

R.  M.  Hochstrosser  and  A.  P.  Marchetti ,  J.  Chem.  Phys.  50,  1727  (1969) 
L.  E.  Harris  H.  S.  Mar  a,  and  S.  P.  McGlynn,  Czech.  J.  Phys.  B20, 

1007  (1970). 
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Via  PUBLICATIONS  PLANNED  AND  IN  PREPARATION 


1.  Vertical  Spectra  for  C02~  and  Electron  Affinity  of  CO2.  E.  Sachs. 


2.  Application  of  the  Differential  Correlation  Energy  Method  to  the 
Calculation  of  the  Electron  Affinity  of  CN  and  C02.  A.  C.  Wahl, 

E.  Sachs,  and  J.  Detrich. 
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8.  The  Differential  Correlation  Energy  Method  for  the  Calculation  of  Electron 
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Int.  J.  Quantum  Chem.  _li  S,  0000  (1980). 
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National  Resource  for  Computation  in  Chemistry,  July,  1980. 

11.  General  Quadratically  Convergent  MCSCF  Theory  in  Terms  of  Reduced  Matrix  Elements 
C.  C.  J.  Roothaan  and  J.  H.  Detrich 
Reviews  of  Modern  Physics 
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VII  INTERACTIONS  AND  ATTENDANCE  AT  MEETINGS 


October  1978  -  J.  Detrich,  Trip  to  Washington  to  discuss  NBS  contract  and 
Wisconsin  contract 

November  28  -  December  1,  1979  -  J.  Detrich,  APS  Meeting,  Division  of  Electron 
and  Atomic  Physics 

May  3-5,  1979  -  J.  Detrich,  Midwest  Theoretical  Chemistry  Conference 
(Poster  Presentation) 

July  11  -  15,  1979  -  J.  Detrich,  OSU  Spectroscopy  Conference  34th  Symposium 
on  Molecular  Spectroscopy  (2  talks) 

August  13  -  18,  1979  -  J.  Detrich,  Gordon  Conference  on  Few  Body  Interactions 

Summer  1979  -  J.  Detrich,  Numerous  discussions  of  new  MCSCF  method  with 
scientists  at  NBS,  NRCC,  and  the  Gordon  Conference 

September  17  &  18,  1979  -  J.  Detrich,  DOD  Physics  Review  -  NAS  Washington, 

D.C. 


March  9  -  15,  1980  -  J.  Detrich,  Sanibel  Sumposium  -  Palmcoast,  Florida 

May  10  &  11,  1980  -  J.  Detrich,  Midwest  Theoretical  Chemistry  Conference, 
Minneapolis  Minnesota 

July  14  -  16,  1980  -  J.  Detrich,  NRCC  -  MCSCF  Workshop  -  College  Station, 
Texas 


May  3-5,  1979  -  A.  C.  Wahl,  Midwest  Theoretical  Chemistry  Conference, 
Chairman  Session  I 

June  1,  1979  -  A.  C.  Wahl,  Discussion  of  current  and  future  status  of  AFOSR 
contract  with  R.  Kelley,  Physics  Directorate  AFOSR,  Washington,  D.C. 

September  17  &  18,  1979  -  A.  C.  Wahl,  DOD  Physics  Review,  NAS  Washington, 

D.C. 

November  5  &  6,  1979  -  A.  C.  Wahl,  DOD  Chemistry  Review,  NAS  Washington, 

D.C. 

March  9  -  15,  1980  A.  C.  Wahl,  Sanibel  Symposium  -  Palmcoast,  Florida 

May  10  &  11,  1980  -  A.  C.  Wahl,  Midwest  Theoretical  Chemistry  Conference, 
Minneapolis,  Minnesota 

June  20  -  25,  1980  -  A.  C.  Wahl,  Spectroscopy  Conference  -  Columbus,  Ohio 

July  14  -  16,  1980  -  A.  C.  Wahl,  NRCC  -  MCSCF  Workshop  -  College  Station, 
Texas 

July  21  -  25,  1980  -  A.  C.  Wahl,  Gordon  Conference  -  Wolfeboro,  New  Hamphsire 
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May  3-5,  1979  -  E.  Sachs,  Midwest  Theoretical  Chemistry  Conference, 

La  Fayette,  Indiana 

July  25  -  27,  1979  -  E.  Sachs,  Workshop  on  Software  Standards  and  Computional 
Chemistry,  Salt  Lake  City,  Utah. 

March  13  -  18,  1978  -  0.  Hopper,  Sanibel  Symporium  -  Palmcoast,  Florida 

March  18  -  20,  1978  -  0.  Hopper,  Eleventh  Midwest  Theoretical  Chemistry 
Conference,  Columbus,  Ohio 

May  28  -  June  2,  1978  -  0.  Hopper,  Twenty-Sixth  Annual  Conference  on  Mass 
Spectrometry  and  Allied  Topics,  St.  Louis,  Missouri 

June  12  -  16,  1978  -  0.  Hopper,  Thirty-Third  Annual  Symposium  on  Molecular 
Spectrometry,  Columbus,  Ohio 

March  11  -  17,  1979  -  D.  Hopper,  Sanibel  Symporium,  Palmcoast,  Florida 

June  3-8,  1979  -  D.  Hopper,  Twenty-Seventh  Annual  Conference  on  Mass 
Spectrometry  and  Allied  Topics. 

July  11  -  15,  1979  -  D.  Hopper,  Thirty-Fourth  Symporium  on  Molecular 
Spectroscopy,  Columbus,  Ohio. 

Sept.  17  -  18,  1979  -  0.  Hopper,  D0D  Physics  Review,  Washington,  D.C. 

Oct.  29  -  Nov.  3,  1979  -  D.  Hopper,  Third  International  Conference  of 
Quantum  Chemistry,  Kyoto,  Japan. 


Over  the  past  year,  A.  C.  Wahl  has  interacted  by  telephone  and  letter  with 
a  number  of  atmospheric  scientists  whose  experimental  work  relates  to  this 
AF0SR  contract  (NBS,  N0AA,  JILA,  SRI). 
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VIII  PROJECT  PERSONNEL 


Below  the  names  and  durations  of  association  with  project  of  personnel  is 
given: 

J.  Detrich  joined  the  project  in  September  1978  and  is  presently  associated 
with  the  project  as-  project  scientist. 

D.  G.  Hopper  joined  the  project  in  March  1978.  His  association  with  the 
project  ended  in  November  1979. 

Ed  Sachs  joined  the  project  in  June  1978  and  left  to  accept  a  permanent 
position  at  Bell  Labs  in  September  1979. 

Frank  Tobin  joined  the  project  in  October  1978,  left  in  June  of  1979  and 
accepted  a  position  at  Johns  Hopkins  University. 

Chris  Wahl  brought  the  project  to  SAI  from  Argonne  oeginning  March  of  1978,  was 
on  a  leave  of  absence  due  to  illness  July  1978  through  June  1979,  and  is  currently 
principal  investigator. 
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IX  Concluding  Remarks 

We  believe  that  this  comprehensive  systematic  project  has  enhanced  our 
knowledge  and  understanding  of  the  atmospheric  triatomic  systems  under 
study.  In  many  cases,  particularly  for  excited  states  and  ions,  new  information 
has  been  generated  which  was  not  and  continues  not  to  be  available  experimentally. 
In  addition,  during  the  course  of  this  work  several  powerful  methods  for 
investigating  molecular  structure  were  developed. 

It  is  our  hope  that  the  Air  Force  will  continue  to  support  this  exciting  work 
which  can  form  the  basis  for  improved  modelling  and  for  new  insights  into 
our  understanding  of  the  chemistry  and  physics  of  the  earth's  atmosphere. 
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APPENDIX  I 


PUBLICATIONS  RESULTING  FROM 
THE  ATMOSPHERIC  PROJECT 


These  publications  report  atmospheric  related  research  performed  with  the 
support  of  the  Air  Force  Office  of  Scientific  Research,  the  Advanced  Research 

T 

Projects  Agency,  the  Aerospace  Research  Laboratories,  the  Atomic  Energy 
Commission,  the  Energy  Research  and  Development  Administration  and  the 
Department  of  Energy  over  the  period  1972  -  1980. 

The  work  represented  by  these  papers  involved  a  number  of  individual  contracts 

involving  A.  C.  Wahl  and  his  collaborators,  ?.  Benioff, 

G.  Das,  J.  Detrich,  W.  England,  P.  Fortune,  D.  Hopper,  R.  Hosteny,  B.  Rosenberg, 
E.  Sachs,  W.  Stevens,  and  C.  W.  Wilson.. 
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TECHNIQUE  PAPERS 

Application  of  Perturbation  Theory  in  Large  Configuration  Interaction 
Calculations 

P.  J.  Fortune  and  B.  J.  Rosenberg 
Chem.  Phys.  Lett.  37,  110-114  (1977) 

A  Modified  Pseopotential  Approach  to  the  Heavy-Atomic  Molecular  Systems: 
Application  to  the  X2I^  +  ,  A2II^,  and  the  A 3 1 1 3^2  States  of  the  HgH  Molecule 

G.  Das  and  A.  C.  Wahl 

J.  Chem.  Phys.  64,  4672-4679  (1976) 

BISON-MC:  A  FORTRAN  Computing  System  for  MCSCF  Calculation  on  Atoms, 

Diatoms  and  Polyatoms.  This  polyatomic  MCSCF-C1  capability 
was  developed  out  of  this  project.  This  has  been  widely 
used  in  polyatomic  MCSCF  calculations  by  many  researchers. 

G.  Das  and  A.  C.  Wahl,  ANL-7955  (1972) 

An  Improved  MCSCF  Method 

C.  C.  J.  Roothaan,  J.  Detrich,  and  D.  G.  Hopper,  I.J.Q.C.S13,93  (1979): 
C.  C.  J.  Roothaan  and  J,  Detrich,  To  be  published. 

.The  Differential  Correlation  Energy  Method  for  Electron  Affinities,  Ionization 
Potentials  and  Vertical  Spectra 

A.  C.  Wahl,  E.  Sachs,  and  J.  H.  Detrich,  International  Journal 
of  Quantum  Chemistry  S14,  0000  (1980). 

Theoretical  Studies  of  Atmospheric  Triatomic  Molecules 

0.  H.  Detrich,  E.  Sachs,  A.  C.  Wahl,  International  Journal 
of  Quantum  Chemistry  S14,  0000  (1980). 
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A  case  History  in  Computer  Resource  Sharing:  Ab  Initio  Calculations  Via  a 
Remote  Control 

O.  G.  Hopper,  P.  J.  Fortune,  A.  C.  Wahl  and  T.  0.  Tiernan 
COMPUTER  NETWORKING  AND  CHEMISTRY,  ACS  SYMPOSIUM  SERIES,  Number  19, 

P.  Lykos,  Editor  (American  Chemical  Society,  Washington,  1975), 
pp  153-182 

The  Multiconfiguration  Self-Consistent  Field  Method 
A.  C.  Wahl  and  G.  Das 

METHODS  IN  ELECTRONIC  STRUCTURE  THEORY,  H.  F.  Schaefer  III,  Editor 
(Plenum,  New  York,  1977 )^,^_p^p  51-73 
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MCSCF  Calculations  on  the  Lowest  Triplet  State  of  H20 
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Chem  Phys.  Lett.  23,  9-12  (1973) 

Theoretical  Studies  of  Atmospheric  Molecules:  SCF  and  Correlated  Potential 
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P.  J.  Fortune,  B.  J.  Rosenberg,  and  A.  C.  Wahl 
J.  Chem.  Phys.  65,  2201-2205  (1976) 
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Recent  applications  of  the  Multiconfiguration  Self-Consistent  Field  Method 
to  Polarizabilities,  Excited  States,  Van  Der  Waals  Forces,  and  Triatomic 
Surfaces 


A.  Karo,  M.  Krauss,  and  A.  C.  Wahl 
Int.  0.  Quantum  Chem.  S7,  143-159  (1973) 

Atmospheric  Studies:  Theoretical  Study  of  the  Ozone  System  (03) 

A.  C.  Wahl  and  C.  W.  Wilson,  Jr. 

Technical  Report  on  ARPA  Order  No.  2022.  (.1972) 

Vibrational  Spectrum  and  Transition  Intensities  of  Ozone 

A.  C.  Wahl,  P.  Benioff,  G.  Das,  and  C.  W.  Wilson,  Jr. 
Technical  Report  on  AFOSR- ISSA-74-004,  (1974) 

CARBON  DIOXIDE  AND  ITS  IONS 


Ab  Initio  Vertical  Spectra  and  Linear  Bent  Correlation  Diagrams  for  the 
Valence  States  of  C02  and  its  Singly  Charged  Ions 
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W.  8.  England,  B.  J.  Rosenberg,  P.  0.  Fortune,  and  A.  C.  Wahl 
J.  Chem.  Phys.  65;  684-691  (1976) 

Theoretical  Studies  of  Atmospheric  Triatomic  Molecules:  Accurate  SCF 
Vertical  Spectrum  for  Valence,  Mixed  Character,  and  Rydberg  States  of 
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W.  B.  England,  W.  C.  Ermler,  and  A.  C.  Wahl 
J.  Chen.  Phys.  6_6_,  2336-2343  (1977) 

Theoretical  Studies  of  Atmospheric  Triatomic  Molecules:  . Ab  Initio  Equations 
of-Motion  Excitation  Energies  for  Valence  States  of  the  Configuration 
lm3g  2— 1 u  in  C02 

W.  B.  England,  D.  Yeager,  and  A.  C.  Wahl 
J.  Chem.  Phys.  66,  2344-2345  (1977) 

Theoretical  Studies  of  Atmospheric  Triatomic  Molecules:  New  Ab  Ini  to 
Results  for  the  -*■  lAu  Vertical  Ordering  in  C02 

W.  B.  England  and  W.  C.  Ermler 
J.  Chem.  Phys.  70,  1711-1719  (1979) 

NI0TR0GEN  DIOXIDE  AND  ITS  IONS 


The  Electronic  Structure  of  Nitrogen  Dioxide.  I.  Multiconfiguration 
Self-Consisten-Field  Calculation  of  the  Low-Lying  Electronic  States 

G.  D.  Gillispie,  A.  U.  Khan,  A.  C.  Wahl,  R.  P.  Hosteny,  and 
M.  Krauss 

J.  Chem.  Phys.  63,  3425-3444  (1975) 

Theoretical  Studies  of  Atmospheric  Molecules:  SCF  and  Correlated  Energy 
Levels  for  the  N02,  N02+  AND  N02"  Systems 

P.  J.  Fortune,  B.  J.  Rosenberg,  W.  B.  England,  and  A.  C.  Wahl 
Theoret.  Chem.  Acta  46,  185  (1977) 
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Ab  Initio  Calculations  of  the  Minimum  Energy  Path  in  the  Doublet  Surface 
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J.  A.  Benioff,  G.  Das,  and  A.  C.  Wahl 
J.  Chem.  Phys.  67,  2449-2462  (1977) 

Ab  Initio  Calculations  of  the  Vertical  Electronic  Spectra  of  N02,  N02+  and 
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P.  A.  Benioff 
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Ab  Initio  Calculations  on  N02  and  N02”:  Optimization  of  Diffuse  Gaussian 
Exponents 

P.  A.  Benioff 

Theoret.  Chem.  Acta.  48,  337  (1978) 

NITROUS  OXIDE  ANO  ITS  IONS 
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D.  G.  Hopper 
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